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Exact  Solution  For  The  Electromagnetic  Fields 
of  a  Uniform  Line  Current  Parallel 
to  a  Flat  Homogeneous  Earth 


1.  INTROD  CTION 

In  his  pioneer  paper,  Carson*  obtained  approximate  solutions  for  essentially 
the  same  problem  addressed  here,  but  the  magnitude  of  the  errors  caused  by  sim¬ 
plifying  assumptions,  such  as  neglecting  displacement  currents,  is  unclear, 
especially  in  cases  of  poor  conductivity  and  high  frequencies.  The  alternate  ap¬ 
proach  adopted  here  is  to  develop  the  exact  solutions  which,  as  might  be  expected, 
are  more  cumbersome.  They  do,  however,  offer  the  possibility  of  numerical 
evaluation  to  a  desired  degree  of  accuracy  in  specific  cases. 

2.  APPROACH 

Figure  1  represents  a  rectangular  coordinate  system  XYZ,  witl  the  XGZ 
plane  coincident  with  the  flat  earth-surface.  A  line- current  of  I  rms  amperes 
flows  in  the  positive  Z-direction  at  height  h,  and  passes  through  the  point  x  -  0, 
y  =  h,  z  =  0.  The  material  of  the  uniform  earth  has  a  dielectric  constant  e  F/m 
and  a  conductivity  of  cr  mho/m.  The  region  y  >  0  is  taken  to  be  free  space  of 

(Received  for  publication  5  April  1982) 

1.  Carson,  J.R.  (1926)  Wave  propagation  in  overhead  wires  with  ground  return, 
Bell  System  Technical  Journal  5:539-554. 
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Figure  1.  Coordinates  in  a  Plane  Transverse  to 
the  Z-axis,  Which  is  Directed  Toward  the  Reader 


dielectric  constant  and  both  media  have  the  magnetic  permeability  of  free 
space,  MKS  units  are  used  throughout. 

The  field  at  an  arbitrary  point  P(x,  y,  z)  above  the  earth  is  regarded  as  consist¬ 
ing  of  three  parts:  (1)  the  fields  which  the  current  would  produce  if  it  were  in  free 
space,  the  earth  being  absent;  (2)  the  fields  of  a  perfect  image  of  the  current  at 
x  =  0,  y  ■  -h,  as  if  the  e?rth  were  perfectly  conducting;  and  (3)  ’’supplementary" 
fields  which  account  for  the  earth  not  being  a  perfect  conductor.  The  fields  inside 
the  earth  have  no  source -singularity  and  are  described  as  just  one  field- complex. 
The  free-space  field  is  discussed  in  Section  3,  the  effects  of  the  image  are  obtained 
in  Section  4,  and  the  supplementary  fields  are  developed  in  Section  5. 


3.  FREE  SPACE  FIELDS 

In  discussing  the  free-space  fields  it  is  convenient  to  use  an  auxiliary  cylindri¬ 
cal  coordinate  system  z)  whose  z-axis  coincides  with  the  line-current,  as 

shown  in  Figure  1.  The  current  I  and  all  its  fields  vary  with  time  t  according  to 
the  factor  e  (which  is  not  explicitly  written  out  in  the  expressions  to  follow) 
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where  w  is  the  angular  frequency  corresponding  to  the  current  frequency  f.  In 
free  space.  Maxwell's  equations  for  the  electric  field  E  and  the  magnetic  field  H 


VXE  =  iupQH 


V  X  H  =  -i  u)  E 


Because  of  symmetry  and  the  absence  of  free  charges,  only  the  field  components 
Ez  and  exist,  and  these  do  not  vary  with  z.  In  cylindrical  coordinates  the 
Maxwell  equations  become 


dT  =  -iu)^oH* 


J~  (r  HJ  =  -iue  E 
r.  otZ  -  v  o  z 

The  Hankel  functions  and  are  related  by  the  recurrence  relations2 

o  1 


dH(1)(u)  m 

-t—  - 


Using  these  it  is  readily  shown  that  Eqs*  (3)  and  (4)  are  satisfied  by 


“4-  «,<»  *rJ 


h*  =  ¥  Hi(1)(kr-) 


where 


k  -  u  /e  U 
V  o  o 


2.  Abramowitz,  A.  ,  and  Stegun,  I.  A.  (1964)  Handbook  of  Mathematical  Functions, 
National  Bureau  of  Standards,  Applied  Mathematics  Series  55,  U.  S.  Govern  - 
ment  Printing  Office,  Washington,  D.  C. 
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where  c  is  the  velocity  of  electromagnetic  waves  in  free  space,  and  \  is  the 
wavelength. 

Other  properties  of  the  Hankel  and  Bessel  functions  of  the  first  kind  (J),  and 
the  second  kind^  (Y)  which  will  be  needed  are: 


H(l)  =  J  +  iY 
o  o  o 


Hl(1>  •  V‘yl 


and  for  u  —0, 


Ho(1>(u) 


,  2i  .  2 

1  -T-  ln  TU 


2i  io  u 

"2  '  Tu  +  T  i  T 


where  y  -  1.  78107. 
Also,  for  v  —  oo. 


H  — 

o 


2  i  (u  -  7r /  4 ) 

1T\)  e 


H  (1)<o>  - 


2  i  (u  -  3?r/4) 

7TV)  e 


Close  to  the  line  current,  kr_  —  0  and 


E  -  '  1  —  In  -1 _ 

z  *4  I1  '  rr  ln  7kr. 


U  _  * 

<t>  2irr  ^ 


The  constants  in  Eqs.  (7)  and  (8)  were  pre-chosen  so  that  Eq.  (12)  agrees  with 
Ampere’s  Law, 


Some  writers  represent  this  function  by  N. 


At  large  distances  from  the  line  current,  kr.  -  »  and 


E 

z 


%ul 

^/Srkr. 


ei(kr.  -  */4) 


(13) 


H 


ei  (kr,  -  ir/4) 


(14) 


Since  the  time  factor  is  e~*u*  and  k  is  a  positive  real  number  these  fields  have  the 
form  of  outgoing  waves,  as  they  should.  Thus  the  fields  given  by  Eqs.  (7)  and  (8) 
meet  all  the  requirements  for  an  exact  solution  for  the  free-space  case.  Also,  for 
kr.  oo. 


(15) 


which  is  the  impedance  relation  for  plane  waves  in  free-space.  The  (radiated) 
power  flowing  through  aim  section  of  a  large  cylinder  is  then 

col2 

2»r.  |Ej  •  |h#|  *  .  (16) 

This  power  loss  is  as  if  the  line -current  encountered  a  resistivity  u  Q  ^  4  S!/m,  and 
it  may  be  seen  from  Eq.  (n)  that  the  in-phase,  or  real  part  of  Ez  at  the  line  cur¬ 
rent  is  negative  and  opposes  the  current  flow.  The  power  required  to  drive  the 
current  against  this  field  is  n0u,'I  /4  W/m  and  is  equal  to  the  radiated  power  The 
imaginary  part  of  E^  has  a  logarithmic  singularity  at  r_  =  0,  as  can  be  seen  from 
Eq.  (11).  This  is  a  consequence  of  having  a  current  of  infinitely  small  cross  sec¬ 
tion;  if  the  current  were  distributed  over  a  finite  cross  section  the  field  would 
remain  finite. 


4.  FIELDS  OVER  A  PERFECTLY  CONDUCTING  EARTH 

The  fields  inside  a  perfectly  conducting  earth  are  of  course  zero,  as  is  the 
tangential  component  of  the  electric  field  at  the  surface.  The  boundary  conditions 
are  met  by  placing  a  180  degree  out -of -phase  image  at  x  =  0,  y  =  -h.  The  fields 
at  a  point  P  above  the  surface  are  then  a  superposition  of  fields  of  the  type  discussed 
in  Section  3.  Thus 


9 


r 


i 


EZ,«  -  [«01)<kr->-H0(1)<kr+)]  (17) 

where  the  subscript  oo  Is  used  to  denote  the  case  of  the  perfectly  conducting  earth, 
and 


r 


± 


+  (y±h)2 


(18) 


Since  the  magnetic  fields  are  in  the  and  directions,  as  shown  in  Figure  1 
and  since 


cos  _ 


h-y 


sin 


x 

r  _ 


and 


(19) 


COS  0-L 


h  +  y 


a  _  X 

K  - 


resolving  the  magnetic  fields  into  x-  and  y-components  gives 

H  =  Kirill  H,(1)(kr  )  +  5L 12.  H ,(1)(kr+)l 
X,  00  4  L  r"  1  "  F+  1  J 

rr  _  ikl  T  x  H  (: 

Hy. 00  '  Hl 


(l)(kr.)  -  ~  H  (1)(kr+) 
r+  1  T  J 


Just  above  the  earth,  as  y-’O,  E  and  H  vanish  while 


(20) 


(21) 


H 


(1) 


H 


_  ikhl  “1 _ \k 


/2~72 

tvx  +  h 


X,  oo 


y  -  o 


(22) 


From  Eq.  (17)  the  electric  field  at  the  current  itself  (x  =  0,  y  =  h)  is 


E  (0,  h)  = 

Z ,  00 


i  ui  r  / 
4-  I  H 

4  L  0 


(1)(0)  -  Ho(1)(2kh) 


(23) 
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RESISTIVITY  RELATIVE  TO 


2 

I 


FREE  SPACE  VALUE 


HEIGHT  OF  WIRE  IN  WAVELENGTHS  -  h/X 

Figure  2.  Radiation  Resistivity  of  a  Wire  Carrying  a  Uniform  Current  as  a  Function 
of  Height  Above  a  Perfectly  Conducting  Flat  Earth 


t2  7t  t2  7t/2 

CJ  I  ry  M  -.U  I  ^ 

— ^ —  /  sinw  (kh  sin  o)  do  =  —  /  (  *  *  cos  #))  do 


V  u>r 

-V“  [  l  -  J  (2kh)J 


where  use  is  made  of  a  standard  integral  form.  This  is,  of  course,  the  same  as 
the  power  expended  in  driving  the  current,  as  discussed  above. 


5,  SUPPLEMENTARY  HELDS  AND  FORMAL  SOLUTION 


Inside  the  earth  the  Maxwell -Ampere  equation  is 

V  X  H  =  (a  -  iue)E  (30) 


instead  of  Eq.  (2),  while  Eq.  (l)  remains  unchanged.  Symmetry  considerations 
show  that  the  only  rectangular  components  of  the  fields  are  £z,  H^,  and  and 
these  do  not  change  with  z.  The  Maxwell  equations  then  reduce  to 


.  dE 
l  _ z 

u*io  9y 


3E 

1  _ z 

uli  Ox 
o 


OH 


OH 

x 

^7 


(cr  -  i  <j  e)  E„ 


(31) 


(32) 


(33) 


It  may  be  shown  by  direct  substitution  that  these  equations  are  satisfied  by  the 
following  expressions: 


00  /2  71 

r  yv i 5  "  ka 

Ez(x.  y)  =  J  G(s)  e 


cos  xs  ds 


Hx(x.y>  = 


•  UJT  f  -  kg2  G(s)  e 


V«2  -  kc 


cos  xs  ds 


H  (x,y)  = 

y  J 


f 

J 


s  G(s)  e 


y\A2  ~K 


sin  xs  ds 


where 


2  2 


=  CJ  €  M  +  1  CT  u  /U 
o  o 


(34) 


(35) 


(36) 


(37) 
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and  G(s)  is  an  arbitrary  function  of  the  variable  of  integration,  to  be  chosen  later 
to  satisfy  the  boundary  conditions.  In  accordance  with  symmetry  requirements 
E  and  H  are  even  functions  of  x,  while  H  is  an  odd  function.  If  a  >  0  the 

z  x  g  ^  y 

quantity  s  -  k  "  is  always  in  the  lower  half  of  the  complex  plane,  and  hence 
/  2  2  ^ 

^Js  -  kg  is  either  in  the  second  or  fourth  quadrant  depending  on  the  choice  of 
sign.  For  Eqs.  (34)*  (35),  and  (36)  to  represent  physically  realistic  fields,  the 
real  part  of  must  be  positive  so  that  the  fields  vanish  as  y  — *  -  oe.  Thus 

-  kg  is  taken  to  be  in  the  fourth  quadrant. 


At  points  above  the  ground,  Eqs. 
satisfied  by  the  fields 


(31),  (32),  and  (33)  with  £  c  £0  anc*  a  ~  0  are 


r  -  /  2"  "T? 

E^Cx,  y)  =  J  A(s)  e  cos  xs  ds  +  E^  { 


(38) 


Hx(x,y)  = 


/  2  2 

— —  f  v42  -  k2  A(s)  e  ^  ~  ^  cos  xs  ds  +  H 

up  j  * 


(39) 


Hy(x,y)  = 


l 

"uu 


/ 


s  A(s)  e 


-yv42  - 1 


sin  xs  ds  +  H 


(40) 


where  E  ,  H  ,  and  H  _  are  the  fields  for  the  case  of  perfect  ground  conduc- 
Z,  00  X,  00  y,  oo  v  b 

tivity  already  discussed  in  Section  4.  The  parts  of  the  fields  containing  the  arbi¬ 
trary  function  A(s)  individually  satisfy  the  Maxwell  equations  if  k  *  u  „  /e_  p  as 
2  2  /2  2  v  o  o 

before.  When  s  >  k  ,  the  positive  real  value  for  vs  -k  must  be  chosen  so 

2  2 

that  the  fields  decrease  with  increasing  y.  If  s  <  k  ,  the  negative  value  for 
/~2  2 

Vs  -k  must  be  chosen  to  have  the  solution  correspond  to  outgoing  waves  rather 
than  physically  unrealistic  incoming  waves.  As  a  reminder  of  this  sign  choice,  an 
asterisk  (  )  is  added  to  the  square  root  sign  in  subsequent  expressions. 

At  the  earth’s  surface  (y  =  0)  the  tangential  components  of  the  electric  and 
magnetic  fields  must  be  continuous.  Continuity  of  requires 


00 

l'z(xt0)  =  J*  A  (s)  cos  xs  ds  - 
o 


/ 


G(s)  cos  >s  ds 


(4  1) 


since  K 


.  x  vanishes  at  y  =  0.  Evidently  A  and  G  are  essentially  Fourier  trans¬ 
forms  of  the  same  even  function  E  (x,0),  so  that 

z 
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A(s)  *  G(s) 


Similarly,  continuity  of  requires 


—i—  f  Vs-k2  A(s)  cosxs  ds  +  H 


x,  oc  UjU 


42  -  kg2  G(s)  cosxs  ds 


On  using  Eqs.  (22)  and  (42)  this  condition  becomes 


^2 - 2  *  /2 - 2  1  khl  H^Mk^+h2 

s  -k  +  V»  -  k  1  A(s)  cosxs  ds  =  - * — ;  2 ~ - 

*  J  vx  +h 


■  i/™’ 


s)  cosxs  ds 


where  by  the  Fourier  integral  theorem  the  (dimensionless)  function 


,s)  “  k  f  - - -  H  J  ^  (  k  Vx2  +Ti2  )  cos  xs  dx 

J  fTy  9  I  \  t 


o  Vx  +  h 


It  follows  that 


u  u;I 

A(s)  =  -  -2— 

7T 


F(h,  s) _ 

/  2  ,2  .  *  / 2  .2 

's  -kg  +  Vs  'k 
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On  consolidating  the  results,  the  fields  above  the  earth  are: 


Ez<*.y> 


•vra 


F(h,  s)  cos  xs 


ds 


‘tT  [Ho(1)(krJ  -  Ho1>(kr+)] 


H  (x.  y)  =  -  —  [  e'y  —  ~  k  V^-k2  F(h,  s)  cos  xs  dg 

x  T  J  +  v^-T 


V®2  -  kg2  +  Vs* 

Hi(1>(kr-)  +  ^Hi(1)(kr+)] 


Hy(x'y)  =  T  / 


-V  *  /2~U 

e  y  _ sF(h,  s)  cos  xs 


/  2  ,  2  .  *  /2~2 
Vs  *  kg  +  Vs 


ds 


X  f  Hi(1)(kr.)  -  £  H1(l)(kr+) 


Below  the  surface  of  the  earth  the  fields  are: 


Ez<x.y>  = 


/2  .  2 

F(h,s)  cos  xs 


^/s2  -  kg2  +  Vs2  -  k2 


ds 


Hx(x,y) 


■t/ 


00  /2  .  2 
Vs  _k 


9  /2  ST 

^/s  -kg  F(h,s)  cos  xs 


ds 
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r 

H  (x,y)  =  H  /  ® _ 8  F(h,s)  sin  xs  dg 

y  ’ 1 

O  * 

On  converting  to  the  dimensionless  parameters: 

u  =  s/k 
x^  =  kx 

yx  =  ky 

hx  =  kh 

Eq.  (45)  becomes 


Ffhj.u)  =  J  -}==. .  '  Hl(1)(^i2  +  hi2  )  cosuxjdxj 

°  VX1  +hl 


and,  for  instance ,  with  y  >  0, 


00  £ 

-y 


/j  i 

e 

— 


'  /2  7 

vu  -  1 


F  (h u)  cos  ux  j 


u2-k,V+^ 


du 


^  U3  I 


V-  [Ho<1)(krJ  *  Ho(1)(kr+>] 


where 


kgW 


E/£o  +  ia/uEo 


(52) 


(53) 


(54) 


(55) 


(5G) 


6.  PLAUSIBLE  APPROXIMATIONS 

Before  proceeding  with  a  numerical  investigation  of  a  particular  case  of  the 
exact  solution,  it  may  be  of  interest  to  mention,  without  justification,  some 


1 

i 


( 


\  \ 
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plausible  approximations.  If  the  real  part  of  the  Hankel  function  H ^  is  ignored, 
and  its  small -argument  approximation 

H,(l)(u)  a  -  (57) 

I 

(see  Eq.  (10)]  is  used  in  the  integral  for  F(h,u)  even  though  the  range  of  integra¬ 
tion  extends  to  infinity,  then 


h,  cos  ux 


F(h  ^  u)  = 


2i  f  nl  06 

vJ 

o  1 


~T~  ^1 


(58) 


This  integral  is  a  standard  form,  giving 


Ffrijiu)  =  -ie 


-hjU 


(59) 


Next  if  W-  1  is  replaced  by  u,  Eq.  (55)  gives 


oo  -(h  +y.)u 

ifi  ul  re  cos  ux,  muI  r  m  n)  -i 

km  »  -V-  /  7s=Tt — 1 *■  -  -J-k"'""-1  -  H„“M 

a/u  -k_/k  +u 

(60) 


J  / 2  .  2/2  . 

o  u  -  k_  /k  +  u 


Finally,  if  the  Hankel  functions  IlJ1^  are  replaced  by  small  argument  approxima¬ 
tions: 


Ho1)(kr±)  a  ‘--T  (n^ 


(see  Eq.  (10))  ,  the  result  is 


iM  uj I 


Ez(x,y)  a  ~2T 


(hj  +  yjhi 


(  /  a/u2  -  kg2  /  k2  +  u 


c.  V*,2  +  (yi+hi)' 

-i-  du  +  in  z  - - — 

x/x2  +  (yt  -h^ 


(61) 
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3 

This  is  similar  to  a  form  used  by  Wait  .or  highly  conducting  earths  and  low 
frequencies.  ^ 


7.  NUMERICAL  EXAMPLE 


Using  the  exact  solution  (Eq.  (55)] ,  the  real  part  of  the  electric  field  at  the 
lin*>  current  itself  (x  =  0;  y  s  h)  will  be  calculated  numerically  to  within  a  few  per¬ 
cent  for  the  following  case: 

f  =  3.  5  X  104  Hz 

A  =  8571.42857  m 

k  =  7.  33C 38287  X  lo'4  m'1 

a  =  10  mho/m 

e/eo  =  10 

a/ut  =  514.2857143 
o 

h  =  10  m 

y j  =  hj  =  7.  33038287  X  10‘3 

kr+  =  2kh  =  0.014660764 

kr  _  =  0 

Then 

Re  I  H^Vkr.)  -  =  JQ(0)  -  JQ(0.  014660764) 

=  5. 38  X  10‘5 

4 

by  interpolation  of  tables,  and 


Re(Ez] 


V"1 


/|i  ul 

(Qr  +  iQ^  F(h  ,u)du  -  (5.  38  X  Iff  D) 


(62) 


^Wait’s  formulation  evidently  contains  a  typographical  error  in  which  his  sym¬ 
bols  x  and  y  were  interchanged  in  three  places. 

3.  Wait,  J.R.  (1961)  On  the  impedance  of  a  long  wire  suspended  over  the  ground, 

Proc.  IRE,  49,  October,  p  1576. 

4.  Cambi,  E.  (1948)  Eleven  and  Fifteen-Place  Tables  of  Bessel  Functions  of  the 

First  Kind  to  all  Significant  Orders,  Dover,  New  York 
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where 


Q  =  Qr  +  iQi 


(63) 


The  functions  and  in  the  interval  0  <  u  <  400  are  illustrated  in  Figure  3. 
Then  if  Fr  and  are,  respectively,  the  real  and  imaginary  parts  of  Ffti^.u), 


Ke{E„ 


M  ul 
O 


w  uu 

f  Q  F  du  -  -  f  Q.  F. 
J  r  r  if  J  ii 


du  +  1.  08  X  10 


-4 


(64) 


Now 


F 

r 


hi  / 


jj(p> 


cos  ux^  dx^ 


(65) 


with 


P 


(66) 


Figure  4  shows  a  graph  of  J  j(p)/p  for  0  <  x1  <  20.  For  large  values  of  xj»  J j(p)/p 
approaches  y/lj  n  cos  (Xj  -  3tt/4)/x^^^  so  that 


f  Ji(p) 

J  Q 

20 


cos  uXj  dxj 


00 


vTo 


0.  36 


for  all  values  of  u.  The  portion  of  the  integral  from  0  to  2  0  is  not  larger  than  about 

.3 

0.  5,  so  it  is  concluded  that  Fr  is  less  than  something  of  the  order  of  or  7  X  10 

for  all  values  of  u,  and  will  be  much  less  than  h^  at  large  values  of  u.  Inspection 

of  Figure  3  shows  that  the  area  under  the  curve  of  Q  (u)  is  about  1,  so  that  the 

r  *  3 

first  term  inside  the  large  bracket  in  Eq.  (64)  is  not  larger  than  about  7  X  10 


20 


2h. 


D(u)  = 


/ 


A  cos  uXj  dx^ 


*  *  +  p 


Evaluating  the  integral  of  Eq.  (68)  as  in  Section  6, 


(69) 


(70) 


Figure  5.  A  Portion  of  the  Func  tion  A(x^) 


to  the  function  D,  for  all  values  of  u.  For  values  of  x1  greater  than  10,  the 
limiting  form 


A  (x  j y/ijr 


_  sin  (x,  -  3 it/ 4)  , 

3r  _ 1 _  4.  1 


represents  A  with  an  error  less  than  about  10  .  The  contribution  to  D  from 


Xj  >  10  is  then  less  than  about 


for  all  values  of  u.  Thus  for  0  5  u  5  100,  L)(u)  ^  2h^/ ir  *  0.  005,  or  less  than 
1  percent  of  e  ^lU  which  is  between  1  and  0.480  in  this  range.  At  larger  values 
of  u,  the  entire  factor  F.  becomes  quite  unimportant  because  the  multiplying  fac¬ 
tor  Q.(u)  in  Eq.  (84)  is  then  much  less  than  1  percent  of  its  initial  value,  as  shown 
in  Figure  3. 

The  product  function  -Qi  F.  in  Eq.  (64)  is  shown  in  Figure  6  for  the  important 
part  of  its  range.  A  numerical  integration  gives 


0.  33 


(74) 


00 

J  -Qi  (u)  F^u)  du  £ 

o 

and  hence 


RelEJ 


V"1 


[i 


(0.  33)  - 


U>"3) 


+  1.  1  X  10 


1 


ii  (J  I 

*  -°-84-%- 


The  corresponding  resistivity  of  the  line  is  0.  84  *iQu>/8 


0.  02  9  ft/m.t 


(75) 


^Based  on  the  previously  mentioned  approximate  formulation  (Section  6), 
Wait**®  gave  two  terms  of  a  series  for  the  resistivity: 


Y  ["  h*  •■•] 


0.200,  and 


For  the  above  numerical  example,  the  quantity  i  Qu  h 

the  corresponding  line  resistivity  is  0.  80  m0u)/8  ft/m. 

5,  Wait,  J.R.  (1974)  Comments  on  'The  horizontal  wire  antenna  over  a  conducting 
or  dielectric  half  space:  current  and  admittance,  1  Radio  Science  9:1165. 
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5  veillance  of,  ground  and  aero*pace  object*,  intelligence  data  t 
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h  compatibility. 
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